Non-Markovian Decay and Lasing Condition in an Optical Microcavity Coupled to a 
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The decay dynamics of the classical electromagnetic field in a leaky optical resonator supporting 
a single mode coupled to a structured continuum of modes (reservoir) is theoretically investigated, 
and the issue of threshold condition for lasing in presence of an inverted medium is comprehensively 
addressed. Specific analytical results are given for a single-mode microcavity resonantly coupled 
to a coupled resonator optical waveguide (CROW), which supports a band of continuous modes 
acting as decay channels. For weak coupling, the usual exponential Weisskopf-Wigner (Markovian) 
decay of the field in the bare resonator is found, and the threshold for lasing increases linearly 
with the coupling strength. As the coupling between the microcavity and the structured reservoir 
increases, the field decay in the passive cavity shows non exponential features, and correspondingly 
the threshold for lasing ceases to increase, reaching a maximum and then starting to decrease as the 
coupling strength is further increased. A singular behavior for the "laser phase transition", which 
is a clear signature of strong non-Markovian dynamics, is found at critical values of the coupling 
between the microcavity and the reservoir. 

PACS numbers: 42.55.Ah, 42.60.Da, 42.55.Sa, 42.55.Tv 
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I. INTRODUCTION. 

It is well known that the modes of an open optical 
cavity are always leaky due to energy escape to the out- 
side. Mode leakage can be generally viewed as due to 
the coupling of the discrete cavity modes with a broad 
spectrum of modes of the " universe" that acts as a reser- 
voir [3, H, Q . From this perspective the problem of es- 
cape of a classical electromagnetic field from an open 
resonator is analogous to the rather general problem of 
the decay of a discrete state coupled to a broad contin- 
uum, as originally studied by Fano ^4,] and encountered 
in different physical contexts (see, e.g., Q). The sim- 
plest and much used way to account for mode coupling 
with the outside is to eliminate the reservoir degrees of 
freedom by the introduction of quasi normal modes with 
complex eigenfrequencies (see, e.g., [H, Q)) in such a way 
that energy escape to the outside is simply accounted for 
by the cavity decay rate 7 (the imaginary part of the 
eigenvalue) or, equivalently, by the cavity quality factor 
Q. This irreversible exponential decay of the mode into 
the continuum corresponds to the well-known Weisskopf- 
Wigner decay and relies on the so-called Markovian ap- 
proximation (see, e.g., Q) that assumes an instantaneous 
reservoir response (i.e. no memory): coupling with the 
reservoir is dealt as a Markovian process and the evo- 
lution of the field in the cavity depends solely on the 
present state and not on any previous state of the reser- 
voir. For the whole system (cavity plus outside), in the 
Markovian approximation the cavity quasi-mode with a 
complex frequency corresponds to a resonance state with 
a Lorentzian lineshape. If now the field in the cavity ex- 
periences gain due to coupling with an inverted atomic 
medium, the condition for lasing is simply obtained when 
gain due to lasing atoms cancels cavity losses, i.e. for 
g ~ 7, where g is the modal gain coefficient per unit 



time [H . More generally, treating the field classically and 
assuming that the cavity supports a single mode, an ini- 
tial field amplitude in the cavity will exponentially decay, 
remain stationary (delta-function lineshape) or exponen- 
tially grow (in the early stage of lasing) depending on 
whether 5 < 7, 5 = 7 or 17 > 7, respectively. In addition, 
since the cavity decay rate 7 increases as the coupling of 
the cavity with the outside increases, the threshold for 
laser oscillation increases as the coupling strength of the 
resonator with the modes of the "universe" is increased. 
It is remarkable that this simple and widely acknowl- 
edged dynamical behavior of basic laser theory, found in 
any elementary laser textbook (see, e.g., [6|), relies on the 
Markovian assumption for the cold cavity decay dynam- 
ics [7]. However, it is known that in many problems deal- 
ing with the decay of a discrete state coupled to a "struc- 
tured" reservoir, such as in photoionization in the vicin- 
ity of an autoionizing resonance 8] , spontaneous emission 
and laser-driven atom d yna mics in waveguides and pho- 
tonic crystals [1, 0, [HI, M E El [H, El [13 > and elec- 
tron transport in semiconductor superlattices [Tsj, the 
Markovian approximation may become invalid, and the 
precise structure of the reservoir (continuum) should be 
properly considered. Non-Markovian effects may become 
of major relevance in pre sence of threshold [l^l or sin- 
gularities [13, E2j EIj Tisl . Ell in the density of states or 
more generally when the coupling strength from the ini- 
tial discrete state to the continuum becomes as large as 
the width of the continuum density of state distribution 
. Typical features of non-Markovian dynamics found in 
the above-mentioned contexts are non-exponential decay, 
fractional decay and population trapping, atom-photon 
bound states, damped Rabi oscillations, etc. Though 
the role of structured reservoirs on basic quantum elec- 
trodynamics and quantum optics phenomena beyond the 
Markovian approximation has received a great attention 
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(see, e.g., Ref.p^ for a rather recent review), at a clas- 
sical level previous works have mainly considered the 
limit of Markovian dynamics developing a formal- 
ism based on quasi-normal mode analysis of the open 
system Q. In fact, in a typical laser resonator made 
e.g. of two-mirrors with one partially transmitting mir- 
ror coupled to the outside open space, the Weisskopf- 
Wigner decay law for the bare cavity field is an excel- 
lent approximation |l| and therefore non-Markovian ef- 
fects are fully negligible. However, the advent of micro- 
and nano-photonic structures, notably photonic crystals 
(PCs), has enabled the design and realization of high- 
Q passive microcavities [20I. l2ll. [2^ [23l [2^ and lasers 
[Hi m, m, 113, ill which can be suitably coupled to 
the outside by means of engineered waveguide structures 
[2l|,[2i,[33,[33,[3l,[3l|. By e.g. modifying some units cells 
within a PC, one can create defects that support local- 
ized high-Q modes or propagating waveguide modes. If 
we couple localized defect modes with waveguides, many 
interesting photon transport effects may occur (see, e.g., 
[291 . [sol . [341 ). Coupling between optical waveguides and 
high-Q resonators in different geometries has been inves- 
tigated in great detail using numerical methods, coupled- 
mode equations, and scattering matrix techniques in the 
framework of a rather general Fano- Anderson-like Hamil- 
tonian [2^, [sO, Ull, [H, [H, [11] . Another kind of light cou- 
pling and transport that has received an increasing atten- 
tion in recent years is based on coupled resonator optical 
waveguide (CROW) structures [H, [13, [M [11 , in which 
photons hop from one evanescent defect mode of a cavity 
to the neighboring one due to overlapping between the 
tightly confined modes at each defect site. The possibility 
of artificially control the coupling of a microcavity with 
the "universe" may then invalidate the usual Markovian 
approximation for the (classical) electromagnetic field de- 
cay. In such a situation, for the passive cavity one should 
expect to observe non-Markovian features in the dynam- 
ics of the decaying field, such as non-exponential decay, 
damped Rabi oscillations, and quenched decay for strong 
couplings. More interesting, for an active (i.e. with gain) 
microcavity the usual condition (7 = 7 of gain/loss bal- 
ance for laser oscillation becomes meaningless owing to 
the impossibility of precisely define a cavity decay rate 7. 
Therefore the determination of the lasing condition for a 
microcavity coupled to a structured reservoir requires a 
detailed account of the mode structure of the universe 
and may show unusual features. 

It is the aim of this work to provide some general insights 
into the classical-field decay dynamics and lasing condi- 
tion of an optical microcavity coupled to a structured 
reservoir, in which the usual Markovian approximation 
of treating the cavity decay becomes inadequate. Some 
general results are provided for a generic Hamiltonian 
model describing the coupling of a single-mode micro- 
cavity with a continuous band of modes, and the effects 
of non-Markovian dynamics on lasing condition are dis- 
cussed. As an illustrative example, the case of a micro- 
cavity resonantly coupled to a CROW is considered, for 



which analytical results may be given in a closed form. 
The paper is organized as follows. In Sec. II a simple 
model describing the classical field dynamics in an ac- 
tive single-mode microcavity coupled to a band of con- 
tinuous modes is presented, and the Markovian dynam- 
ics attained in the weak coupling regime is briefly re- 
viewed. Section III deals with the exact dynamics, be- 
yond the Markovian approximation, for both the passive 
(i.e. without gain) and active microcavity. In partic- 
ular, the general relation expressing threshold for laser 
oscillation is derived, and its dependence on the cou- 
pling strength between the microcavity and the reservoir 
is discussed. The general results of Sec. Ill are special- 
ized in Sec. IV for the case of a single-mode microcavity 
tunneling-coupled to a CROW, and some unusual dy- 
namical effects (such as "uncertainty" of laser threshold, 
non-exponential onset of lasing instability and transient 
non-normal amplification) are shown to occur at certain 
critical couplings. 



II. MICROCAVITY COUPLED TO A 
STRUCTURED RESERVOIR: DESCRIPTION OF 
THE MODEL AND MARKOVIAN DYNAMICS 



A. The model 



The starting point of our anal ysis is provided by a 
rather general Hamiltonian model [29l [soj describing the 
interaction of a localized mode \a) of a resonator system 
(e.g. a microcavity in a PC) with a set of continuous 
modes \ujfj_) of neighboring waveguides with which the 
resonator is tunneling-coupled. We assume that the mi- 
crocavity supports a single and high-Q localized mode 
of frequency uja, and indicate by 7^ and g the intrin- 
sic losses and gain coefficients of the mode. The intrin- 
sic losses 7i account for both internal (e.g. absorption) 
losses and damping of the cavity mode due to coupling 
with a "Markovian" reservoir (i.e. coupling with modes 
of the universe other than the neighboring waveguides). 
The modal gain parameter g may be provided by an in- 
verted atomic or semiconductor medium hosted in the 
microcavity. Since we will consider the microcavity op- 
erating below or at the onset of threshold for lasing, as in 
Refs. [13, [1E| the modal gain parameter g is assumed to be 
a constant and externally controllable parameter; above 
threshold an additional rate equation for g would be obvi- 
ously needed depending on the specific gain medium (see, 
for instance, [40| ) . Dissipation and gain of the microcav- 
ity mode are simply included in the model by adding 
a non-Hermitian term Hnh to the Hermitian part of 
the Hamiltonian. The full Hamiltonian H then reads 
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H = Ho + Hint + Hnh: where % 



Ho = u;a\a){a\+^ duj^uj^\ujf,)(ujf,\, 
Hint = J dLjf^[Kf_,{ujfj,)\ujfj,){a\ + h.c 



Hnh = i(5-7i)|a)(a|, 



(la) 

(lb) 
(Ic) 



with {a\a) = 1, = '5^y5(w^ - w^), (a|w^) = 0, 

and h = 1. The coefficients k^{lu^) describe the direct 
couphng between the locahzed mode \a) of the microcav- 



ity and the propagating modes \uj^) in the continuum, 
whereas A is a dimensionless parameter that measures 
the strength of interaction (A ^ for a vanishing inter- 
action). If we write the state {ip) as 



^ Ca{t)\a) / dujf,Cf,{uj^,t)\ 



(2) 



the following coupled-mode equations for the coefficients 
Ca(t) and c^{u!fi, t) are readily obtained from the equation 
id\^p)/dt = H\iP): 



iCa{t) = {ljJa+ ig ~ ili)Ca{t) + j dcj^K* (tJp)Cp (cJ^ , i) , 



I 



(3a) 
(3b) 



where the dot stands for the derivative with respect 
to time t. Note that the power of the microcavity 
mode is given by |ca(t)P, whereas the total power of 
the field (cavity plus structured reservoir) is given by 

P{t) = |ca(t)P + E;./rf^/.|cMK,*)l'- The threshold 
condition for lasing is obtained when an initial pertur- 
bation in the system does not decay with time. From 
Eqs.(l3a| and (j3bp the following power-balance equation 
can be derived 

dP 

^ = (5-70K|^ (4) 

from which we see that |cap for any g < 7^, so that 
the threshold g = gth for laser oscillation satisfies the 
condition gth > 7i, as expected. 

B. Weak coupling limit: Markovian dynamics 

The temporal evolution of the microcavity-mode am- 
plitude Ca{t) and the condition for laser oscillation can 
be rigorously obtained by solving the coupled-mode equa- 
tions pap and (j3b[) by means of a Laplace transform anal- 
ysis, which will be done in the next section. Here we 
show that, in the weak coupling regime (A — s- 0) and for 
a broad band of continuous modes, coupling of the cav- 
ity mode with the neighboring waveguides leads to the 
usual Weisskopf-Wigner (exponential) decay. Though 
this is a rather standard result (see, e.g. Q) and ear- 
lier derived for a standard Fabry-Perot laser resonator 
in Ref. [l| using a Fano diagonalization technique, for the 
sake of completeness it is briefly reviewed here within 
the model described in Sec. II. A. If the system is initially 
prepared in state |a), i.e. if at initial time t = there is 
no field in the neighboring waveguides and Ca(0) 7^ 0, an 



integro-differential equation describing the temporal evo- 
lution of cavity mode amplitude Ca {t) at successive times 
can be derived after elimination of the reservoir degrees 
of freedom. A formal integration of Eas.(|3bp with initial 
condition 0^(0;^, 0) = yields 

c^{ujf_,,t) ^ -iXKf^{uj^) / dt'ca{t')exp[-iuj^{t - t')]. 
Jo 

(5) 

After setting Ca{t) — A{t) exp{—iujat), substitution of 
Eq.® into Ea. ([5a|) yields the following exact integro- 
differential equation for the mode amplitude A{t) 

A={g-j,)A- f dTG{T)A{t~T), (6) 
Jo 

where G(t) is the reservoir response (memory) function, 
given by 

G(r) = A^^ / du;^|K^(cj^)pexp[-z(cj^ -t^a)T]. (7) 

Equation ^ clearly shows that the dynamics is not a 
Markovian process since the evolution of the mode am- 
plitude at time t depends on previous states of the reser- 
voir. Nevertheless, if the characteristic memory time 
is short enough (i.e., the spectral coupling coefficients 
broad enough) and the coupling weak enough such that 
|^/yl|r,„ ^ 1, we may replace Eq.® with the following 
approximate equation 

A^{g-^,)A-A{t) f dTG{T) ~ {g-j,)A~ilBMAR)A, 
Jo 

(8) 

where 

(7fl + iAr) = / dTG{T) (9) 
Jo 
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for t Tm- In this limit, the dynamics is therefore 
Markovian and the reservoir is simply accounted for by 
a decay rate and a frequency shift Ar. Using the 
relation 



lim / (irexp(- 



-iuJT) = 7r(5(ti>) — iV 



OJ 



(10) 



from Eq.([7]) the following expressions for the decay rate 
-fa and the frequency shift can be derived 



^ J UJa-i^u 



(11) 

(12) 



The dynamics of the cavity mode field in the Markovian 
approximation is therefore standard: an initial field am- 
plitude in the cavity will exponentially decay, remain sta- 
tionary (delta-function lineshape) or exponentially grow 
(in the early stage of lasing) depending on whether 5 < 7, 
g = 7 or 5 > 7, respectively, where 7 = 7i -f is the to- 
tal cavity decay rate. The threshold for laser oscillation 
is therefore simply given hy gtn — li + 1r, i-e. 



gth 



= 7^-1- ttA^ ^ |k^(cJq)P. 



(13) 



III. FIELD DYNAMICS BEYOND THE 
MARKOVIAN LIMIT: GENERAL ASPECTS 

Let us assume that the system is initially prepared in 
state |a), i.e. that at initial time i = there is no field 
in the neighboring waveguides [c^{uj^,Q) ~ 0] whereas 
Cq(0) = 1. The exact solution for the field amplitude 
Cq {t) of the microcavity mode at successive times can be 
obtained by a Laplace- Fourier transform of Eas. (|3ap and 
(I3bp . Let us indicate by Cq(s) and c^{llj^, s) the Laplace 
transforms of Ca{t) and c^{uj^,t), respectively, i.e. 



Ca{s) 



dt Cait) exp(— si) 



(14) 



and a similar expression for c)i(aj^,s). From the power 
balance equation (|4]), one can easily show that the in- 
tegral on the right hand side in Eq. p^ converges for 
Re(s) > 77, where 77 = for 5 — 7i < or 77 = g — 7^ for 
g — "fi > 0. The field amplitude Ca{t) is then written as 
the inverse Laplace transform 

Ca{t) = -r^ / ds Ca{s) exp(si) (15) 

where the Bromwich path B is a vertical line Re(s) = 
const > 77 in the half-plane of analyticity of the trans- 
form, and Co(s) is readily derived after Laplace transform 
of Eqs. (|3a|) and pb]) and reads 



Ca(s) 



is — LUa — ig' — 5](s) 



(16) 



In Eq. ([T6)) . g' = g — % is the effective gain parameter 
and S(s) is the self-energy function, which is expressed 
in terms of the form factor 



, V(lu) 
dio 

zs — w 



(17) 



where I?(a;) is the reservoir structure function, defined 
by 



(18) 



In writing Eg. pTj) . we assumed that the spectrum of 
modes of the waveguides (to which the microcavity is 
coupled) shows an upper and lower frequency limits ui 
and UJ2- We will also assume that ^{uj) does not show 
gaps, i.e. intervals with V = 0, inside the range (a;i,<x'2). 
The assumption of a finite spectral extension for the con- 
tinuous modes is physically reasonable and is valid for e.g. 
PC waveguides or CROW. In addition, in order to avoid 
the existence of bound states (or polariton modes) for 
the passive microcavity coupled to the structured reser- 
voir, we assume that I?(w) vanishes at the boundary of 
the band, precisely we require that ^{uj) ^ {uj — uj 1^2)^^ 
as Lo —f wi,2, with 61^2 > 0. This condition, which will 
be clarified in Sec.III.A, is a necessary requirement to 
ensure that the field amplitude Ca (t) fully decays toward 
zero for g' = 0. 

The temporal evolution of Ca{t) is largely influenced by 
the analytic properties of Ca{s); in particular the occur- 
rence of a singularity (pole) at s = Spoie with Re(spo/e) > 
may indicate the onset of an instability, i.e. a las- 
ing regime. The self-energy function S(s) [Eq. (jl7p ]. and 
hence Ca{s), are not defined on the segment of the imag- 
inary axis s = —iu! with tui < uj < W2, si^2 = —i^i,2 
being two branch points. In fact, using the relation 

lim = r (-] T inSiuj), (19) 

p_»o+ u}±tp \^ J 

from Eg. pT)) one has 

E(s = -iw±0+) = A(w) =F«7rX>(w), (20) 
(lui < u) < 102), where we have set 



A(w) =V duj' 



(21) 



To further discuss the analytic properties of Ca{s) and 
hence the temporal dynamics of Ca{t), one should distin- 
guish the cases of passive {g' = 0) and active {g' > 0) 
microcavities. 



A. The passive microcavity 

Let us first consider the case of = 0, i.e. of a pas- 
sive microcavity with negligible internal losses. In this 
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A(Q) 



straight line 



(a) 



A(Q) 




(b) 



FIG. 1: Graphical determination of the roots of Eq. (|23p be- 
low (a), and above (b) the critical coupling. In (b) the full 
Hamiltonian H = Hq + Hint has discrete eigenvalues corre- 
sponding to bound modes. 



case the full Hamiltonian is Hermitian {Hmh — 0), and 
therefore the analytic properties of Ca{s) and spectrum 
oi H = Hn + Hint are ruled as follows (see, for instance, 
SHIill,!!!): (i) The eigenvalues lu oi H are real- valued 
and comprise the continuous spectrum uji < lu < uj2 
of unbounded modes and up to two isolated real-valued 
eigenvalues, outside the continuous spectrum from either 
sides, which correspond to possible bound (or polariton) 
modes [l^; (ii) The isolated eigenvalues are the poles 
of Ca(s) on the imaginary axis outside the branch cut 
— W2 < Im(s) < —LOi; (iii) Ca{s) is analytic in the full 
complex plane, apart from the branch cut and the two 
possible poles on the imaginary axis corresponding to 
bound modes; (iv) In the absence of bound modes Ca{t) 
fully decays toward zero, whereas a limited (or fractional) 
decay occurs in the opposite case. 

From Eq. p^ . the poles s = — ifi of Ca(s) outside the 
branch cut are found as solutions of the equation: 



fl — UJ„ 



duj- 



— to 



i.e. [see Eq.dST])]: 



(22) 



(23) 



with the constraint > ^2 or < wi [43]. A graphical 
solution of Eq. as intersection of the curves VL — oja 
and A(ri) is helpful to decide whether there exist poles 
of Ca(s), i.e. bound modes (see Fig.l). To this aim, note 
that A(f]) > and dl^/dU, < for f7 > cja, A(0) < 
and dA/dV, < for fl < uji, and limj2^±oo A(r2) — 0^. 
Therefore, Eq.(P5|) does not have solutions outside the 
interval {uJi,u;2) provided that A (0)2) < ^^2 — and 
A(a;i) > LUi — LOa [Fig. 1(a)]. Such conditions require at 
least that tOa be internal to the band (cL;i,a;2), i.e. that 
the resonance frequency uja of the microcavity be embed- 
ded in the continuum of decay channels, and that 'D(uj) 
vanishes as a power law at the boundary lo = lui and 
u! = UJ2, i-e. that I?(w) ^ {uj — ti'i.2)''^'^ as w — s- uJi^2 
for some positive integers 61 and 62- In fact, if 'D{uj) 



does not vanish as a power law at these boundaries, 
one would have A(f2) — > ±00 as —^ 0^2, tfi. Even 
though r>(tt') vanishes at the boundaries, as the couphng 
strength A is increased either one or both of the condi- 
tions A (0^2) > ^2 — and A(a;i) < uji — uja can be 
satisfied [Fig. 1(b)], leading to the appearance of either 
one or two bound states. The coupling strength at which 
a bound state starts to appear is referred to as criti- 
cal coupling. Below the critical coupling [Fig. 1(a)], for 
the passive microcavity Ca(s) does not have poles and a 
complete decay of Ca{t) is attained. However, owing to 
non-Markovian effects the decay dynamics may greatly 
deviate from the usual Weisskop-Wigner exponential de- 
cay. The exact decay law for Ca{t) is obtained by the 
inverse Laplace transform Ea. (|15p . which can be eval- 
uated by the residue method after suitably closing the 
Bromwich path B with a contour in the Re(s) < half- 
plane (see, e.g. [l pp.220-221, and [ill. El]). Since the 
closure crosses the branch cut —1^2 < Im(s) < —uji on the 
imaginary axis, the contour must necessarily pass into the 
second Riemannian sheet in the section of the half-plane 
with —UJ2 < Im(s) < — wi, whereas it remains in the first 
Riemannian sheet in the other two sections Im(s) > —lui 
and Im(s) < —IV2 of the Re(s) < half-plane. To prop- 
erly close the contour, it is thus necessary to go back and 
turn around the two branch points of the cut at s = —iuji 
and s — —iuj2, f^oUowing the Hankcl paths hi and /i2 as 
shown in Fig. 2. Note that, while Cq(s) is analytic in the 
first Riemannian sheet for Re(s) < 0, the analytic contin- 
uation c^^(s) of Cq(s) from the right [Re(s) > 0] to the left 
[Re(s) < 0] half-plane across the cut has usually a simple 
pole at s = Sp with Re(sp) < and —uj2 < Im(sp) < —lui 
(see Fig.2). Since c^^(s) = i/[is - Ua - S"(s)] with 
S"(s) = I](s) - 2TTiV{is) [see Eq.®], the pole Sp is 
found as a solution of the equation 



isp — LOa — 5](sp) + 2iTi'D(iSp) — 0, 



(24) 



I.e. 



i"fp + Ap — dio — 



Ap - i^p- uj 
+2TriV{LJa + Ap - i-fp) = (25) 



where we have set 



-7p - ILJa - iAp. 



(26) 



After inversion, we then find for Ca (t) the following decay 
law 

Ca{t) =Zexp[-7pi-i(wa + Ap)t]+C(t), (27) 

where Z is the residue of c^^ (s) at the pole Sp, and C{t) is 
the contribution from the contour integration along the 
Hankel paths hi and ft-2 (see Fig.2): 



1 



S — — iuJi 



C(t) = — / ds[ci'is)-Ca{s)]eMst) + 



1 

27ri 



s— 0— iaJ2 



ds [ci\s) - Ca{s)\ exp(si).(28) 



S — — OQ — IUJ-2 
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FIG. 2: Integration contour used to calculate the inverse 
Laplace transform of Ca(s). The bold solid line on the imag- 
inary axis is the branch cut. The integration along the solid 
(dashed) curves is made on the first (second) Riemannian 
sheet of Ca{s). Sp is the pole of Ca{s) on the second Rieman- 
nian sheet in the Re(s) < half-plane. 



The cut contribution C(t) is responsible for the appear- 
ance of non-exponential features in the decay dynamics, 
especially at short and long times; for an extensive and 
detailed analysis we refer the reader to e.g. Refs.[4ll.[43|: 
examples of non-exponential decays will be presented in 
Sec. IV. We just mention here that, in the weak coupling 
limit (T) 0), from Eq. ((25|) one has that 7p and Ap are 
small, and thus using Eq. (fT9|) we can cast Eq. ([25|) in the 
form 



«7p ■ 



, V(lu) 

OLJ 



+ TTiV{uJa) ^ (29) 



from which we recover for the decay rate 7p and frequency 
shift Ap of the resonance the same expressions and 
Afl as given by Eqs. pTj) and in the framework of 
the Weisskopf-Wigner analysis. In the strong coupling 
regime, close to the boundary of appearance of bound 
modes, the decay strongly deviates from an exponential 
law at any time scale, with the appearance of typical 
damped Rabi oscillations (see e.g. Ref. [^], pp. 249- 
255). 



B. Microcavity with gain: lasing condition 

Let us now consider the case of a microcavity with 
gain, i.e. g' > 0. In this case, one (or more) poles Sp 
of Ca{s) on the first Riemannian sheet with Re(s) > 
may appear as the modal gain g' is increased, so that the 
mode amplitude Ca(t) will grow with time, indicating the 
onset of an instability. In this case, the Bromwich path 
B should be closed taking into account the existence of 
one (or more than one) pole in the Re(s) > plane, 
as shown in Fig. 3. For the case of a simple pole Sp = 
— 7p — itUa — iAp, the expression (|27p for the temporal 
evolution of Ca{t) is therefore still valid, where now 7p < 
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FIG. 3: (a) Deformation of the Bromwich path for inverse 
Laplace transformation with one pole Sp on the Re(s) > 
half-plane (unstable state), (b) Corresponding integration 
contour used to calculate the inverse Laplace transform. The 
integration along the solid (dashed) curves is made on the 
first (second) Riemannian sheet of Ca{s). 



and Ap are found as a solution of the equation [compare 
with Eq.jSSI)] 



hp-ig + 



Ap- 



duj- 



V{uj) 



«7p - OJ 



0. (30) 



As a rather general rule, it turns out that, as g' is in- 
creased, the pole Sp of cl^{s), which at 5' = lies in the 
Re(s) < plane, crosses the imaginary axis in the cut 
region. This crossing changes the decay of Ca{t) into a 
non-decaying or growing behavior, and thus it can be as- 
sumed as the threshold for laser oscillation. The modal 
gain at threshold, 5^^, is thus obtained from Eq. (PD|) by 
setting 7p = 0~, i.e. 

- ig't,, + Ap- A{uja + Ap) + inV{uja + Ap) = 0, (31) 
where we used Eq. ([2T|) and the relation 



V 



duj- 



dio- 



Via;) 



Ap — uj 



Ap + iO+ - UJ 

— iTTV{uJa + Ap 



(32) 
(33) 



Therefore the threshold for laser oscillation is given by 

5t/i = 7»+7rI?(w„ + Ap), (34) 

where Ap (the frequency shift of the oscillating mode 
from the microcavity resonance frequency uja) is implic- 
itly defined by the equation 



duj- 



(35) 



i.e. f2osc — = A(r2osc) with ^osc = + Ap. It should 
be noted that, under the conditions stated in Sec. III. A 
ensuring that for the passive microcavity no bound modes 
exist, Eq. (|55|) admits of (at least) one solution for tJa-f Ap 



inside the range (u;i,a;2). The simplest proof thereof can 
be done graphically [see Fig. 1(a)] after observing that 
^2 ~ > A(cj2) and Ui — Ua < A(wi). 
The rather simple Eq. (j34p provides a generalization of 
Eq. (fT3)) for the laser threshold of the active microcavity 
beyond the Markovian approximation and reduces to it 
in the limit Ap ~ 0. The frequency shift Ap, however, 
can not be in general neglected and may strongly affect 
the value of gth in the strong coupling regime. In fact, 
for a small coupling of the microcavity with the struc- 
tured reservoir (A 0), the shift Ap can be neglected 
and therefore gth increases with A according to Eq. p^ . 
However, as A is further increased up to the critical cou- 
pling condition, the shift Ap is no more negligible, and 
the oscillation frequency f2osc = + Ap at lasing thresh- 
old is pushed toward the boundaries oji or u}2, where 'D{ijj) 
and thus g' th vanish. In fact, as A is increased to reach 
the minimum value between A/ // defined by the relation 



^I,II 



{UJI 



UJa 



UJl,2 - OJ 



(36) 



one has fiosc 1^1,2, and hence gth 7i- Therefore, as 
gth initially increases from 7^ as the coupling strength is 
increased from A = 0, it must reach a maximum value 
and then start to decrease until reaching again the 7.; 
value as A approaches the critical value (A/ or A//). As 
the increase of gth with A in the weak coupling regime is 
simply understood as due to the acceleration of the decay 
of the microcavity mode into the neighboring waveguides, 
the successive decreasing of gth is related to the appear- 
ance of a back-coupling of the field from the continuum 
(waveguides) into the microcavity mode, until a bound 
state is formed at the critical coupling strength. 
As a final remark, it should be noted that the precise 
dynamical features and the kind of instability at lasing 
threshold may depend on the specific structure function 
I?(w) of the reservoir. In particular, anomalous dynami- 
cal features may occur at the critical coupling regime, as 
it will be shown in the next section. 



IV. AN EXACTLY-SOLVABLE MODEL: THE 
COUPLING OF A MICROCAVITY WITH A 
COUPLED RESONATOR OPTICAL 
WAVEGUIDE 

To clarify the general results obtained in the previous 
section, we present an illustrative example of exactly- 
solvable model in which a single-mode and high-Q mi- 
crocavity is tunneling-coupled to a CROW structure 
[36, 37, 38, 39], which provides the non-markovian decay 
channel of the microcavity. In a CROW structure, pho- 
tons tunnel from one evanescent defect mode of a cavity 
to the neighboring one due to overlapping between the 
tightly confined modes at each defect site, and therefore 
memory effects are expected to be non-negligible when- 
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FIG. 4: Schematic of a microcavity (M) tunneling-coupled 
to either one (a) or two (b) cavities of a coupled-resonator 
optical waveguide. Plot (c) shows a schematic of a microcav- 
ity coupled with a CROW in the configuration (b) realized 
on a PC planform made of a square lattice of air holes with 
a one-dimensional chain of defects patterned along the lattice 
(Ref.Hg 



ever the coupling rate of the microcavity with the CROW 
becomes comparable with the CROW hopping rate. 



A. The model 

The schematic model of a microcavity tunneling- 
coupled to a CROW is shown in Fig. 4 for two typical 
configurations. The CROW consists of a chain of equally- 
spaced optical waveguides [36, [H, [1^, supporting 
a single band of propagating modes, and the microcav- 
ity is tunneling-coupled to either one [Fig. 4(a)] or two 
[Fig.4(b)] cavities of the CROW. For the sake of defi- 
niteness, we will consider the coupling geometry shown 
in Fig. 4(b), though similar results are obtained for the 
single-coupling configuration of Fig. 4(a). 

The microcavity and the CROW can be realized on a 
same PC planform (see, e.g., [10, E^): the CROW is 
simply obtained by a one-dimensional periodic array of 
defects, placed at distance d and patterned along the lat- 
tice to form resonant cavities with high-Q factors. The 
microcavity is realized by one defect in the array, say 
the one corresponding to index n = 0, which can have 
a resonance frequency uoa different from that of adjacent 
defects and placed at a larger distance do > d than the 
other cavities [see Fig. 4(c)]. The CROW supports a con- 
tinuous band of propagating modes whose dispersion re- 
lation, in the tight-binding approximation, is given by 

M 



uj{k) — luq — 2KC0s(fcc?), 



(37) 
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where k is the hopping amplitude between two consecu- 
tive cavities of the CROW, d is the length of the unit cell 
of the CROW, k is the Bloch wave number, and luq is the 
central frequency of the band. The resonance frequency 
of the microcavity is assumed to be internal to the 
CROW band, i.e. ojq — 2k < ojm < i^o + 2k. The micro- 
cavity is tunneling-coupled to the two adjacent cavities 
of the CROW, and we denote by kq the hopping ampli- 
tude. The ratio kq/k and the position of ujm inside the 
CROW band can be properly controlled by changing the 
geometrical parameters of the defects and the ratio do /d- 
In particular, in the limiting case where the microcavity 
has the same geometry and distance of the other CROW 
cavities, one has kq = k and LUm = i-^o- An excellent 
and simple description of light transport in the system is 
provided by a set of coupled-mode equations for the am- 
plitudes a„ of modes in the cavities (see, e.g., [s^, EH]) 

idn = -K(a„+i + a„_i) (|n| > 2) (38a) 

icL-i — — Ka_2 — kqCq (38b) 

ica = -Ko{a^i + ai) + {uJa + i9)ca (38c) 

idi = —Ka2 — KQCa (38d) 

where Ca is the amplitude of the microcavity mode, g is 
its effective modal gain per unit time, and uja — i-^m — t^o 
is the frequency detuning between the microcavity reso- 
nance frequency ujm and the central frequency luq of the 
CROW band. For e.g. a CROW built in a GaAs-based 
PC with a square lattice of air holes in the design of 
Ref. [i^ , a typical value of the cavity coupling coefhcicnt 
turns out to be k ~ 700 - 800 GHz and luq/k^Sx 10^ 
at the Aq = 850 nm operation wavelength. Note that in 
writing Eqs.(38), we have neglected the internal losses of 
the CROW cavities; a reasonable value of the Q-factor 
for a reahstic microcavity is Q = uJo/{2jioss) ~ 10^ [23 |. 
which would correspond to a cavity loss rate jioss 1 
GHz to be added in Eqs.(38). This loss rate, however, 
is about two-to-three orders of magnitude smaller than 
the cavity coupling coefhcient k, and therefore on a short 
time scale non-Markovian dynamical effects should be ob- 
served even in presence of CROW losses. The effects of 
reservoir (CROW) losses will be briefly discussed at the 
end of the section. 

To study the temporal evolution of an initial field in the 
microcavity, Eqs.(38) are solved with the initial condi- 
tion a„(0) = and Ca(0) = 1. An integral representation 
for the solution of Eqs.(38) might be directly derived in 
the time domain by an extension of the technique de- 
scribed in Refs.[46l. [47j. where a system of coupled- mode 
equations similar to Eqs.(38), but in the conservative (i.e. 
g = 0) case, was considered. However, we prefer here to 
formally place Eqs.(38) into the more general Hamilto- 
nian formalism of Sec.H and then use the Laplace trans- 
form analysis developed in the previous section to obtain 
the temporal evolution for Ca(t). To this aim, in Ap- 
pendix we prove that Ca (t) may be obtained as a solution 
of the following equations, which have the canonical form 



(3) with a simple continuum of modes acting as a decay 
channel 



ica{t) — {uja + ig)ca + X / dLUKf^{uj)c{uj,tj39a) 

ic{uj,t) = LUc{uJ,t) + XKfj_{uj)Ca{t) (39b) 



2k 



with 



Ak^(w) = Ko 



1 



\2kJ 



1/4 



(40) 



Note that the reservoir structure function for this model, 
defined for uji < lu < lu2 with lui = —2k and lu2 = 2k, is 
simply given by 



P(^) 



2kq 



1 - 



(41) 



With this reservoir structure function, the self-energy 
[Eq. fTT]) ] can be calculated in an exact way and reads 



E(.) 



\ K 



2 r 



(42) 



The function A(ijj), as defined by Eq. pT|) . then reads 



{kq/k)'^UJ 



\uj\ < 2k 





-4k2 


'C^2 


-4k2 



Note that the coupling strength between the microcavity 
and the CROW is determined by the ratio Kq/k, the limit 
Kq/ K^ corresponding to the weak coupling regime. 



B. The passive microcavity: from exponential 
decay to damped Rabi oscillations 

Let us consider first the case g = 0. The conditions 
for the non-existence of bound modes, i.e. for a complete 
decay of Ca{t), are cl'2 — > A((jJ2) and wi — Wa < A(a;i) 
(see Sec.HLA), which using Eq. (P5|) read explicitly 



~ 2k ~ \ K 



(44) 



Note that, as a necessary condition, this relation im- 
plies that \uja\ < 2k and (kq/k)'^ < 1. Note also that 
the the critical coupling regime is reached at (kq/k) = 
\/l — |wa|/ (2k). For a coupling strength [kq/k) above 
such a value, the decay of Ca{t) is imperfect due to the 
existence of bound modes between the microcavity and 
the CROW; this case will not be considered here further. 
The temporal decay law for the mode amplitude Ca{t) 
can be generally expressed using the general relation 
((27)) . which highlights the existence of the exponential 
(Wcisskopf-Wigner) decaying term plus its correction due 
to the contribution of the Hankel paths. Perhaps, for the 
microcavity-CROW system it is more suited to make the 
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FIG. 5: Integration contour used for the inverse Laplace 
transform in the passive microcavity-CROW system. 



J 



inverse Laplace transform on the first Riemannian sheet 
of Ca{s) by closing the Bromwich path B with a semi- 
circle with radius i? ^ cx) in the Re(s) < half-plane 
after excluding the branch cut from the domain by the 
contour a as shown in Fig. 5. Since in this case there are 
no singularities of Cq(s), we simply obtain 



Cait) 



1 

2^ 



ds 



exp(si) 



(45) 



which, using Eg . ([20|) . reads explicitly 

exp(— iwt) 



Cait) 



2^ 



duj 



exp(— iwi) 



duj- 



S](-iw- 0+)_ 

'D{ljj) exp(— icijt) 



For the microcavity-CROW model, one then obtains 



(46) 



Ca{t) 



1 



2 /-Zk 



duj- 



exp(— za;^)^/! — (cj/2k)2 



Stt «3 7_2« {{uj/2n) [1 - [k^/kY] - {L0a/2n)Y + (ko/«)*[1 - c^V(4aj2)] 



(47) 



r 



The integral on the right hand side in Eq. (|47|) can be variable uj 
written in a more convenient form with the change of 

I 



-2k cos Q, yielding 



Ca{t) 



dQ 



(fco/«;)^ sin^ Q exp{2iKt cosQ) 



[{ua/^n) + COS Q — (ko/k)2 cos Q] + (kq / k)'^ sin^ I 

I 



(48) 



In this form, the integral can be written [4a | as a series of 
Bessel functions of first kind and of argument 2Kt (Neu- 
mann series). Special cases, for which a simple expression 
for Ca{t) is available, are those corresponding to Wa = 
and Kg = K, for which 



Cait) = M2Kt), 

and to Wa = and kq = n/^/2, for which 



Ca{t) = 



Jl{2Kt) 



Kt 



(49) 



(50) 



Note that the former case corresponds to a critical 
coupling regime, where Ca(s) has two singularities at 
s = ±2iK -I- 0+. The residues of Ca{s) at these singu- 
larities, however, vanish, and therefore the field Ca{t) 
fully decays toward zero with an asymptotic power law 



^ In general, an inspection of the singularities of 

the Ca(s) reveals that, for Wa 7^ 0, at the critical coupling 
strength (kq/k) — \/l — |wa|/(2K) the Laplace trans- 
form Ca(s) has one singularity at either Sp — 2iK + 0^ 
or Sp = —2iK + 0+ of type Ca(s) ~ — Sp. 

The asymptotic decay behavior of Cq (t) at long times can 
be determined by the application of the method of the 
stationary phase to Ea. (|15)) . One then finds that at the 
critical coupling the field Ca{t) decays toward zero with 
an asymptotic power law ~ whereas below the 

critical coupling the decay is faster with an asymptotic 
decay - 1/^3/2. 

Typical examples of non-exponential features in the de- 
cay process as the coupling strength is increased are 
shown in Fig. 6 for Wa = 0. The curves in the figures 
have been obtained by a direct numerical solution of 
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Normalized Time t'=2Kt 

FIG. 6: Decay of the mode amplitude |ca(f)| in a passive 
microcavity-CROW system for coa = and for increasing val- 
ues of coupling strength: (a) kq/h — 0.2, (b) kq/k = 0.707, 
and (c) kq/h = 1 (critical coupling). 

Eqs.(38). Note that, as for weak coupling the exponential 
(Weisskopf-Wigner) decay law is retrieved with a good 
approximation [see Fig. 6(a)], as the coupling strength 
kq/k is increased the decay law strongly deviates from an 
exponential behavior. Note in particular the existence of 
strong oscillations, which are fully analogous to damped 
Rabi oscillations found in the atom-photon interaction 
context 0. For loq ^ 0, the oscillatory behavior of the 
long-time power-law decay is less pronounced and may 
even disappear (see Ref . ^46] ) . 

C. Microcavity with gain 




FIG. 7: Behavior of normalized threshold gain gth/i'2ii) ver- 
sus the coupling strength (kq /i^)^ for a few values of the ratio 

LUa/{2K,). 
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Normalized Time t '=2k( Normalized Gain g/(2K) 



FIG. 8: (a) Behavior of mode amplitude |ca(t')l versus nor- 
malized time t' = 2K.t for {no/f^f = 0.8, u}a/{2K.) = 0.18, 
and for increasing values of normalized gain (//(2k). (b) Be- 
havior of normalized growth rate versus normalized gain for 
{no/nf = 0.8 and uja/{2K) = 0.18. 



Let us consider now the case g > 0. In order to 
determine the threshold for laser oscillation, we have to 
distinguish three cases depending on the value of the 
coupling strength kq/k. 

(i) basing condition below the critical coupling. In 
this case, corresponding to kq/k < \/l — |wa|/ (2k), the 
threshold for laser oscillation is readily obtained from 
Eqs.dSl, jHSl), gl]) and The frequency nose of the 

oscillating mode is given by fiosc = Wa/[1 — {kq/k)'^], 
and the gain for laser oscillation is thus given by 

The typical behavior of normalized threshold gain 
gth/i'^K)'^ versus the coupling strength {kq/k) is shown 
in Fig. 7. Note that, according to the general analysis of 
Sec. III. B, the threshold for laser oscillation first increases 
as the coupling strength is increased, but then it reaches 
a maximum and then decreases toward zero as the crit- 
ical coupling strength is attained. At g = gth, Ca{s) has 



a simple pole at s = Sp = —iQosc + 0^, whereas as g is 
increased above gth the pole Sp invades the Re(s) > 
half-plane. Therefore, the onset of lasing is character- 
ized by an amplitude |ca(i)| which asymptotically decays 
toward zero for g < gth, reaches a steady-state and non- 
vanishing value at g = gth (the field does not decay nor 
grow asymptotically) , whereas it grows exponentially (in 
the early lasing stage) for g > gth with a growth rate 
a{g) — Re(sp) (see Fig. 8). This instability scenario is 
the usual one encountered in the semiclassical theory of 
laser oscillation as a second-order phase transition [48j . 
However, the temporal dynamics at the onset of lasing 
shows unusual oscillations [see Fig. 8(a)] which are a sig- 
nature of non-Markovian dynamics. In addition, as in 
the Markovian limit the growth rate a should increase 
linearly with g — gth, in the strong coupling regime the 
growth rate a shows near threshold an unusual non-linear 
behavior, as shown in Fig.8(b). 

(a) Lasing condition at the critical coupling with 7^ 0. 
A different dynamics occurs when the coupling strength 
k/kq reaches the critical limit kq/k = y^l — \u!a\/(2K). 
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FIG. 9: Same as Fig. 8, but for parameter values {kq/k)^ — 
0.8 and uia/i^n) — 0.2 (critical coupling). Note that in this 
case there exists no lasing threshold in the traditional sense. 



As discussed in Sec. IV. B, at g = the Laplace transform 
Ca(s) has a singularity at either Sp — 2iK or Sp = —2iK, 
however Sp is not a simple pole and Ca (t) asymptotically 
decays toward zero. For uja 0, i.e. for {k/kq) < 1, 
as g is increased just above zero Ca(s) shows a simple 
pole with a growth rate a — Re(sp) > which slowly 
increases with g at the early stage, as shown in Fig. 9. In 
the figure, a typical temporal evolution of Ca{t) is also 
shown. Note that in this case there is not a value of 
g for which the field amplitude Ca(t) does not grow nor 
decay, i.e. the intermediate situation shown in Fig. 8 (a) 
is missed in Fig. 9(a): for g = the amplitude decays, 
however for = 0^ it always grows exponentially. The 
transition describing the passage of laser from below to 
above threshold in the linear stage of the instability is 
therefore quite unusual at the critical coupling. 

(iii) Lasing condition at the critical coupling withoja = 0. 
A somewhat singular behavior occurs at the critical cou- 
pling when LDa — 0, and therefore kq/k = 1. This case 
corresponds to consider a periodic CROW in which one 
of the cavities is pumped and acts as the microcavity in 
our general model. For = and kq/k — 1, the Laplace 
transform Ca(s) is explicitly given by 



Cq(s) 



1 



-g + \/s2 + 4k2 



(52) 



To perform the inversion, one needs to distinguish four 
cases. 

(a) g — 0. For 5 = 0, the field Ca{t) decays according to 



Ca{t) = Jo(2Ki) 



(53) 



as shown in Sec. IV. B. 



(b) < g < 2k. In this case Ca{s) has two simple poles on 
the first Riemannian sheet at si^2 — ±i\/4j<F~~g^ + 0+. 
The inversion can be performed by closing the Bromwich 
path B with the contour shown in Fig. 10, where along 
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FIG. 10: Integration contour used to calculate the inverse 
Laplace transform for uJa ~ 0, kq/k — 1 and for < g/{2K) < 
1. The integration along the solid (dashed) curves is made on 
the first (second) Riemannian sheet of Ca{s). si,2 are the two 
poles of Ca{s) on the imaginary axis inside the cut. 



the dashed curves the integrals are performed on the 
second Riemannian sheet. One then obtains 



Ca{t) = 



25 



: sm 



(V4K2-.g2t) +C{t) (54) 



where the first term on the right hand side in the equa- 
tion arises from the residues at poles si.2, whereas C{t) is 
the contribution from the contour integration along the 
Hankel paths hi and /12, which asymptotically decays to- 
ward zero as t ~f 00. Note that, after an initial transient, 
the amplitude |ca(t)| steadily oscillates in time with fre- 
quency -^/i^^"^^^ and amplitude 2g/ y^An'^ — g^ . Note 
also that the amplitude and period of oscillations diverge 
as the modal gain g approaches 2k~ . 

(c) g — 2k. In this case, Ca{s) has a single pole of second- 
order in s = 0^ , and therefore to perform the inversion it 
is worth separating the singular and non-singular parts 

of Ca{s) as 



Ca{s) 



4k 



(55) 



where /(s) has no singularities on the imaginary axis. 
After inversion one then obtains 



Gait) = 4k< + 



1 

2^ 



duj /{—iuj + 0~^) exp{~iujt), (56) 



where the second term on the right-hand side in the 
above equation asymptotically decays toward zero. 
Therefore, we may conclude that at 5 = 2k the mode 
amplitude Ca{t) is dominated by a secular growing term 
which is not exponential. 

(d) g > 2k. In this case, Ca(s) has an unstable 
simple pole at Sp — {g^ ~ 'in'^y/^, and therefore the 
solution Ca{t) grows exponentially with time. 

The dynamical scenario described above for = 
and kq/k = 1 is illustrated in Fig.ll. Note that in this 
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FIG. 11: Behavior of mode amplitude Ca{t') versus normal- 
ized time t' = 2Mit for u>a = 0, /to/'* = 1 (critical coupling) 
and for increasing values of normalized gain: (a) g/{2K) — 0, 
(b) g/{2K) = 0.2, (c) ff/(2Ac) = 0.95, (d) ff/(2^c) = 1, and (d) 
5/(2^) = l.l. 



case there is some uncertainty in the definition of laser 
threshold, since there exists an entire interval of modal 
gain values, from g = 0+ to 5 = 2k~ , at which an initial 
field in the cavity does not grow nor decay. 

As a final comment, we briefly discuss the effects 
of internal losses of the CROW cavities, which have 
been so far neglected, on the temporal evolution of 
the mode amplitude Ca{t). In the case where all the 
cavities in the CROW have the same loss rate ^loss, 
the temporal evolution of Ca{t) is simply modified by 
the introduction of an additional exponential damping 
factor exp(-7;ossOi i-S- Ca(i) Cait) exp{-jiosst) ■ This 
additional decay term would therefore shift the threshold 
for laser oscillation to higher values and, most impor- 
tantly for our analysis, it might hinder non-Markovian 
dynamical effects discussed so far. However, for a small 
value of ^loss/i^ (e.g. Jioss/k ~ 0.01 for the numerical 
values given in Ref.fi^l), non-Markovian effects should 
be clearly observable in the transient field dynamics for 
times shorter than ~ i/jioss- As an example, Fig. 12 
shows the dynamical evolution of the mode amplitude 
|ca(t)| for the same parameter values of Fig.ll, except 
for the inclusion of a CROW loss rate jioss = 0.01k. 



It is worth commenting on the dynamical behavior of 
Fig. 12(d) corresponding to g — 2k. In this case, using 
Eq. ([55]) and disregarding the decaying term on the right 
hand side in Eq. (|56p . one can write 

Ca{t) 'iKtex.p{~-fiosst)- (57) 

Note that in the early transient stage the initial mode 
amplitude stored in the microcavity linearly grows as in 
Fig. 11(d), however it reaches a maximum and then it fi- 
nally decays owing to the prevalence of the loss-induced 
exponential term over the linear growing term. There- 
fore, though the microcavity is below threshold for oscil- 
lation as an initial field in the cavity asymptotically de- 
cays to zero, before decaying an initial field is subjected 
to a transient amplification. The maximum amplification 
factor in the transient is about 2k/^ioss, and can be 
therefore relatively large in high-Q microcavities. Such a 
transient growth despite the asymptotic stability of the 
zero solution should be related to the circumstance that 
for 7^ the system (38) is non- normal [1^: though its 
eigenvalues have all a negative real part, the system can 
sustain a transient energy growth. The transient am- 
plification shown in Fig. 12(d) is therefore analogous to 
non-normal energy growth encountered in other hydro- 
dynamic [13, [m, [111 and optical [HI, [H, [5^ systems and 
it is an indicator of a major sensitivity of the system to 
noise. 



V. CONCLUSIONS 

In this work it has been analytically studied, within a 
rather general Hamiltonian model [Eqs.(l)], the dynam- 
ics of a classical field in a single-mode optical microcavity 
coupled to a structured continuum of modes (reservoir) 
beyond the usual Weisskopf-Wigner (Markovian) approx- 
imation. Typical non-Markovian effects for the passive 
microcavity are non-exponential decay and damped Rabi 
oscillations (Sec. III. A). In presence of gain, the general 
condition for laser oscillation, that extends the usual 
gain/loss rate balance condition of elementary laser the- 
ory, has been derived (Sec. III. B), and the behavior of 
the laser threshold versus the microcavity-reservoir cou- 
pling has been determined. The general results have been 
specialized for an exactly-solvable model, which can be 
implemented in a photonic crystal with defects: an opti- 
cal microcavity tunneling-coupled to a coupled-resonator 
optical waveguide (Sec. IV). A special attention has been 
devoted to study the transition describing laser oscilla- 
tion at the critical coupling between the cavity and the 
waveguide (Sec. IV. C). Unusual dynamical effects, which 
are a clear signature of a non-Markovian dynamics, have 
been illustrated, including: the existence of a finite in- 
terval of modal gain where the field oscillates without 
decaying nor growing, the gain parameter controlling the 
amplitude and period of the oscillations; a linear (in- 
stead of exponential) growth of the field at the onset of 
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FIG. 12: Same as Fig. 11, but in presence of CROW losses 
i-yicss/n = 0.01). 



imation and the canonical formulation for the decay of 
a discrete state into a continuum provided by Eqs.(39). 
To this aim, let us first note that, owing to the inversion- 
symmetry of the initial condition a_„(0) = a„(0) = 
(n ^ 0), it can be readily shown that the solution 
a„(i) maintains the same symmetry at any time, i.e. 
o--nit) = an{t) for t > 0. Let us then introduce the 
continuous function of the real-valued parameter Q 



<j){Q, t) = ^ an{t) sin(nQ), 



(Al) 



ra=l 



where Q is taken inside the interval [0, tt]. Using the 
relation 



dQ sm{nQ) sin(m(3) = —Sm.n {m, n > 1) (A2) 
2 



the amplitudes a„ of modes in the CROW are related to 
the continuous field <j) by the simple relations 



a„(i) 



2 r 



dQ 4){Q,t)sm{nQ) (A3) 



(n > 1). The equation of motion for (/) is readily obtained 
from Eqs.(38) and reads 

i— = -2k cos{Q)(j) - Ko sin(Q)ca (A4) 



instability for laser oscillation; and the existence of tran- 
sient (non-normal) amplification of the field below laser 
threshold when intrinsic losses of the microcavity are 
considered. It is envisaged that, though non-Markovian 
effects are not relevant in standard laser resonators in 
which the field stored in the cavity is coupled to the 
broad continuum of modes of the external open space 
by a partially-transmitting mirror [l[ , they should be ob- 
servable when dealing with high-Q microcavities coupled 
to waveguides, which act as a structured decay channel 
for the field stored in the microcavity. 



whereas the equation for Ca, taking into account that 
a_i -I- fli = 2ai = (4/7r) dQ(f>{Q, t) sa\{Q), can be cast 
in the form: 

ica{t) = {uj^^+ig)c,{t)-— f dQ 0(g,i)sin(Q). (A5) 



By introducing the frequency uj of the continuum 

w = -2kcos(Q) (A6) 

and after setting 



APPENDIX A 

In this Appendix it is proved the equivalence between 
coupled-mode equations (38) in the tight-binding approx- 



/ 2 1 

C{uj,t) = -\ <l>{uJ,t) T—r. 

V TTK^' '[1-C^2/(2k)2]1/4' 



(A7) 



one finally obtains Eqs. (|39ap and (|39b[) given in the text. 
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